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Abstract 



We show that the problem of finding an e-approximatc Nash equilibrium of an 
n x n two-person games can be reduced to the computation of an (e/n) 2 -approximate 
market equilibrium of a Leontief economy. Together with a recent result of Chen, 
Deng and Teng, this polynomial reduction implies that the Leontief market exchange 
problem does not have a fully polynomial-time approximation scheme, that is, there is no 
algorithm that can compute an e-approximate market equilibrium in time polynomial 
in m, n, and 1/e, unless PPAD C P, We also extend the analysis of our reduction 
to show, unless PPAD C RP, that the smoothed complexity of the Scarf's general 
fixed-point approximation algorithm (when applying to solve the approximate Leontief 
market exchange problem) or of any algorithm for computing an approximate market 
equilibrium of Leontief economies is not polynomial in n and 1/c, under Gaussian or 
uniform perturbations with magnitude a. 
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1 Introduction 



A Leontief economy |231 I28j with m divisible goods (or commodities) and n traders is 
specified by a pair of non-negative m x n matrices 1 E = (ejj) and D = (dij). After 
trading, the exchange of goods can also be expressed by a non-negative m x n matrix 

X = (Xi : j). 

• The matrix E is the endowment matrix of the traders, that the amount of 
commodity i that trader j initially has. We can assume, without loss of generality, 
that there is one unit of each type of good. With this assumption, the sum of every 
row of E is equal to 1. 

• The matrix D is the demand matrix or the utility matrix. It defines n utility functions 
u\, . . . ,u n , one for each trader. For each exchange X of goods, let be the j th column 
of X, then 



mm 

i 



X i,j 
dij 



is the Leontief utility of trader j. 

The initial utility of trader j is Uj(ej) where is the j th column of E. The individ- 
ual objective of each trader is to maximize his or her utility. However, the utilities that 
these traders can achieve depend on the initial endowments, the individual utilities, and 
potentially, the (complex) process that they perform their exchanges. 

In Walras' pure view of economics the individual objectives of traders and their 
initial endowments enable the market to establish a price vector p of the goods in the 
market. Then the whole exchange can be conceptually characterized as: the traders sell 
their endowments - to obtain money or budgets - and individually optimize their objectives 
by buying the bundles of goods that maximize their utilities. 

By selling the initial endowment e^, trader j obtains a budget of (e^lp) amount, where 
(ej|p) denotes the dot-product of these two vectors. The optimal bundle for Uj is a solution 
to the following mathematical program: 

maxtij(xj) subject to (xj|p) < (ej|p) . (1) 

A solution to Equation Q is referred to as an optimal demand of trader j under prices p. 
The price vector p is a Walrasian equilibrium, an Arrow-Debreu equilibrium, or simply an 
equilibrium of the Leontief economy (E, D) if there exists optimal solution X — (xi, • • • , x n ) 
to Equation Q such that 



< 1, (2) 



where 1 is the m-dimensional column vector with all ones. This last constraint states that 
the traders' optimal demands can be met by the market. 

In other words, an equilibrium price vector essentially allows each trader to make indi- 
vidual decision without considering others' utilities nor how they achieve their objectives. 



See Section[5]for the basic notations used in this paper. 
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The computation of a market equilibrium is a fundamental problem in modern economics 
23 as Walrasian equilibria might provide useful information for the prediction of market 
trends, in the decision for future investments, and in the development of economic policies. 
So a central complexity question in Leontief market exchange problem is: 

Question 1 (Polynomial Leontief?). Is the problem of computing an equilibrium of a 
Leontief economy in P? 

So far no polynomial-time algorithm has been found for this problem. In practice, 
one may be willing to relax the condition of equilibria and considers the computation of 
approximate market equilibria. For example, Scarf j^2] developed a general algorithm for 
computing approximate fixed points and equilibria. Recently, Deng, Papadimitriou, Safra 
|B] proposed a notion of an approximate market equilibrium as a price vector that allows 
each trader to independently and approximately optimize her utilities. 

For any e > 0, let OPTj e (p, e^) be the set of e- approximately optimal vectors in 
for Equation ([T|). that is, the set of all Xj satisfying 



Then p is an e-approximate equilibrium 2 of a Leontief economy (E, D) if there exists 
e- approximately optimal solution X = (xi, . . . , x n ) with Xj 6 OPT Jj(: (p, ej) such that 



Question 2 (Fully Polynomial Approximate Leontief?). Can an e-approximate equi- 
librium of a Leontief economy with m goods and n traders be computed in time polynomial 
in m, n, and 1/e? 

The combination of two recent results greatly dashed the hope for a positive answer 
to Question ^ Codenotti, Saberi, Varadarajan, and Ye (Hj gave a polynomial-time reduc- 
tion from two-person games to a special case of the Leontief economy. In a remarkable 
breakthrough, Chen and Deng [3j subsequently proved that the problem of finding a Nash 
equilibrium of a two-person game is PPAD-complete 3 . 

2 One can of course define a stronger notion of approximate equilibria: p is an e-strictly approximate 
equilibrium if there exists X = (xi, . . . , x„) such that 



It is easy to show that every (e/2)-approximate Walrasian equilibrium p is an e-strictly approximate Wal- 
rasian equilibrium, by dividing its associated exchange X by a factor of (1 — e/2). Therefore, in the remainder 
of this paper, we will, without loss of generality, stay with the less restrictive notion of approximation equi- 
libria. 

3 We refer the readers who are not familiar with the complexity class PPAD to the paper by Papadimitriou 
|21| and also to the subsequent papers on the PPAD-completeness of normal games |5] 1210 ED El- This class 
includes several important search problems such as various discrete fixed-point problems and the problem 
of finding a Nash equilibrium of r-person games for any fixed r. 



(xj|p) < (1 + e) (ej\p) , and 



%(x?) > (1 -e)u(Xj-), Vx^- : (x^-|p) < (e^p) 




(3) 



(x^p) < (e f [p), 

Wj(xj) > (1 - e)w(Xj), Vxj- : (xj |p) < ( ej |p) , 



anal 
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In this paper, we show that there is no fully polynomial-time approximation scheme 
for Leontief economies unless PPAD C P. Hence it is unlikely that Question [21 has a 
positive answer. By analyzing the numerical properties of the reduction of Condenotti, 
Saberi, Varadarajan, and Ye [3], we prove that the problem of finding an e-approximate 
Nash equilibrium of an n x n two-person games can be reduced to the computation of 
an (e/n) 2 -approximate equilibrium of a Leontief economy with In goods and 2n traders. 
This polynomial relationship between approximate Nash equilibria of bimatrix games and 
approximate Walrasian equilibria of Leontief economies is significant because it enables us 
to apply the recent result of Chen, Deng, and Teng to show that finding an approximate 
market equilibrium with only 0(log n)-bits of precision is as hard as finding an exact market 
equilibrium, which in turn is as hard as finding a Nash equilibrium of a two-person game 
or a discrete Brouwer fixed point in the most general settings. 

We also consider the smoothed complexity of the Leontief market exchange problem. In 
the smoothed model introduced by Spielman and Teng [21] , an algorithm receives and solves 
a perturbed instances. The smoothed complexity of an algorithm is the maximum over its 
inputs of the expected running time of the algorithm under slight random perturbations 
of that input. The smoothed complexity is then measured as a function of both the input 
length and the magnitude a of the perturbations. An algorithm has smoothed polynomial- 
time complexity if its smoothed measure is polynomial in n, the problem size, and ml/a 



In the smoothed model for Leontief economies, we start with a pair ofmxn matrices 



The smoothed complexity of the Leontief exchange problem (E, D) is then measured by the 
expected complexity of finding an equilibrium of the Leontief economy (E, D). 

The following has been an open question in the smoothed analysis of algorithms. 

Question 3 (Smoothed Polynomial Leontief?). Can an equilibrium of a Leontief econ- 
omy be computed in smoothed time polynomial in m, n, and 1/a? 

Can an e- equilibrium of a Leontief economy be computed in smoothed time polynomial 
in m, n, 1/e and 1/a? 

A concrete open question has been whether the smoothed complexity of the classic 
Scarf's general fixed-point approximation algorithm [22] is polynomial for solving the Leon- 
tief market exchange problem. 

By refining our analysis of the reduction from the two-person games to Leontief economies, 
we show it is unlikely that Scarf's algorithm has polynomial smoothed complexity for com- 
puting an approximate equilibrium of Leontief economies. In particular, we prove that, 
unless PPAD C RP, the problem of finding an (approximate) equilibrium of a Leontief 
economy is not in smoothed polynomial time. 

2 Notations 

We will use bold lower-case Roman letters such as x, a, hj to denote vectors. Whenever a 
vector, say a £ M. n is present, its components will be denoted by lower-case Roman letters 



[21] EE] 
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with subscripts, such as a±, . . . , a n . Matrices are denoted by bold upper-case Roman letters 
such as A and scalars are usually denoted by lower-case Roman letters. The (i,j) th entry 
of a matrix A is denoted by ctij. We use aj to denote the i th column of A. 
We now enumerate some other notations that are used in this paper. 

• R+: the set of m-dimensional vectors with non-negative real entries. 

• Rj^™: the set of all m x n matrices with real entries between a and b. For example, 
Rp.gi" is the set of non- negative matrices with entries at most 2. 

• P ra : the set of all vectors x in n dimensions such that Y27=l Xi = ^ an< ^ Xi — ^ ^ or an 
1 < % < n. 

• (a|b): the dot-product of two vectors in the same dimension. 

• ||x|| p : the p-norm of vector x, that is, |^f|) 1//p and ((x^ = max, \xi\. 

3 Approximate Nash Equilibria of Two-Person Games 

The non- zero-sum two-person game or the bimatrix game is a non-cooperative game between 
two players |18l 1151 116j , the row player and the column player. If the row player has m pure 
strategies and the column player has n pure strategies, then their payoffs are given by a 
pair ofmxn matrices (A, B). 

A mixed row strategy is a vector x G W m and a mixed column strategy is a vector 
y G P n . The expected payoffs to these two players are respectively x T Ay and x T By. A 
Nash equilibrium is then a pair of vectors (x* G F m ,y* G P n ) such that for all pairs of 
vectors x G P m and y G P", 

(x*) T Ay* > x T Ay* and (x*) T By* > (x*) T By. 

Every two-person game has at least one Nash equilibrium ^Hj. But in a recent breakthrough, 
Chen and Deng [3j proved that the problem of computing a Nash equilibrium of a two-person 
game is PPAD-complete. 

One can relax the condition of Nash equilibria and considers approximate Nash equilib- 
ria. There are two possible notions of approximation. 

Definition 3.1 (Approximate Nash equilibria). An e-approximate Nash equilibrium 

of game (A, B) is a pair of mixed strategies (x*, y*) ; such that for all x, y G P n , 

(x*) T Ay* > x T Ay* - e and (x*) T By* > (x*) T By - e. 

Definition 3.2 (Relatively Approximate Nash equilibria). An e-relatively-approximate 
Nash equilibrium of game (A, B) is a pair of mixed strategies (x*,y*), such that for all 

xjeP", 

(x*) T Ay* > (1 - e)x T Ay* and (x*) T By* > (1 - e)(x*) T By. 
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Note that the Nash equilibria and the relatively-approximate Nash equilibria of a two- 
person game (A, B) are invariant under positive scalings, i.e., the bimatrix game (c\ A, C2B) 
has the same set of Nash equilibria and relatively-approximate Nash equilibria as the bi- 
matrix game (A,B), as long as c\,C2 > 0. However, each e-approximate Nash equilibrium 
(x, y) of (A, B) becomes a c- e-approximate Nash equilibrium of the bimatrix game (cA, cB) 
for c > 0. 

On the other hand, Nash equilibria and approximate Nash equilibria are invariant under 
shifting, that is, for any constants c\ and C2, the bimatrix game [c\ + A,C2 + B) has the 
same set of Nash equilibria and approximate Nash equilibria as the bimatrix game (A, B). 
However, shifting may not preserve the relatively-approximate Nash equilibria. 

Thus, we often normalize the matrices A and B so that all their entries are between 
and 1, or between -1 and 1, in order to study the complexity of approximate Nash equilibria 

[miii. 

Recently, Chen, Deng, and Teng [I] proved the following result. 

Theorem 3.3 (Chen-Deng- Teng). The problem of computing an 1/ra 6 - approximate Nash 
equilibrium of a normalized n x n two-person game is PPAD- complete. 

As pointed out in [I], the 6 in the exponent of the above theorem can be replaced by 
any positive constant. One can easily derive the following corollary. 

Corollary 3.4 (Relative approximation is also hard). It remains PP AD -complete to 
compute a l/n e ^ -relatively-approximate Nash equilibrium of an n x n two-person game. 

4 Leontief Market Equilibria: Approximation and Smoothed 
Complexity 

In this and the next sections, we analyze a reduction tt that transforms a two-person game 
(A, B) into a Leontief economy (E,D) = 7r(A,B) such that from each (e/n) 2 -approximate 
Walrasian equilibrium of (E, D) we can construct an e-relatively- approximate Nash equi- 
librium of (A, B). 

We will also consider the smoothed complexity of Leontief economies. To establish 
a hardness result for computing an (approximate) market equilibrium in the smoothed 
model, we will examine the relationship of Walrasian equilibria of a perturbed instance 
(E, D) of (E, D) and approximate Nash equilibria of (A, B). In particular, we show that if 
the magnitude of the perturbation is <r, then we can construct an (e + n 1 ' 5 y / o r )-relatively- 
approximate Nash equilibrium of (A, B) from each (e/re) 2 -approximate equilibria of (E, D). 

Because the analysis needed for approximate Walrasian equilibria is a special case (with 
a = 0) of the analysis for the smoothed model, we will write only one proof for this general 
case, and present it in the next section. 

In this section, we discuss the reduction from two-person games to Leontief economies, 
connect the smoothed model with approximation, and present the main theorems of this 
paper. 

4.1 Approximate Market Equilibria of Leontief Economy 

We first introduce a form of approximate market equilibria that is easier for the analysis of 
the reduction between game equilibria and market equilibria. 



6 



Let D be the demand matrix and E be the endowment matrix of a Leontief economy 
with m goods and n traders. Given a price vector p, trader j can obtain a budget of 
(edp) by selling the endowment. By a simple variational argument, one can show that the 
optimal demand Xj with budget (ej|p) satisfies Xij/ckj = xyjldyj for all i and i' with 
dij > and &y j > 0. Thus, under the price vector p, the maximum utility that trader j 
can achieve is if (ej|p) = 0, and (e,-|p) / (dj|p) otherwise. Moreover, in the latter case, 

= di,j (( e j|p) / (dj|p)). Let u = (u±, . . . , u n ) denote the vector of utilities of the traders. 
Then p is a Walrasian equilibrium price if 

P>0, u; = -^-r, and Du < 1. (4) 
(di|p) 

In the remainder of this paper, we will refer to a pair of vectors (u, p) that satisfies 
Equation @ as an equilibrium of the Leontief economy (E,D). Then, an e-approximate 
equilibrium of the Leontief economy (E, D) is a pair of utility and price vectors (u, p) 
satisfying: 

Ui > (1 — e) (ej|p) / (dj|p) , Vi. - All traders are approximately satisfied. 
Ui < (1 + e) (e$|p) / (dj|p) ,Vi. - Budget constraints approximately hold. 
Du < (1 + e) • 1. - The demands approximately meet the supply. 

Note that if (u, p) is an equilibrium of (E, D), so is (u, ap) for every a > 0. Similarly, if 
(u, p) is an e-equilibrium of (E, D), so is (u, ap) for every a > 0. Thus, we can normalize p 
so that HpI!} = 1. In addition, for approximate equilibria, we assume without loss of gener- 
ality that u and p are strictly positive to avoid division-by-zero since a small perturbation 
of an approximate equilibrium is still a good approximate equilibrium. 

4.2 Reduction from NASH to LEONTIEF 

Let (A,B) be a two-person game in which each player has n strategies. Below we assume 
A G K^ x 2 ™ and B £ K 5 ?. We use the reduction introduced by Codenotti, Saberi, Varadara- 
jan, and Ye [H] to map a bimatrix game to a Leontief economy. This reduction constructs a 
Leontief economy with (E,D) = 7r(A, B) where the endowment matrix is simply E = l2 n , 
the (2n) x (2n) identity matrix and the utility matrix is given by 




(E,D) is a special form of Leontief exchange economies [281 15] . It has 2n goods and 2n 
traders. The j th trader comes to the market with one unit of the j'^-good. In addition, the 
traders are divided into two groups Ai = {1, 2, n} and M = {n + 1, 2n}. Traders in 
A4 only interests in the goods associated with traders in M and vice versa. 

Codenotti et al j5] prove that there is a one-to-one correspondence between Nash equi- 
libria of the two person game (A, B) and market equilibria of Leontief economy (E, D). It 
thus follows from the theorem of Nash |19| llSj. that the Leontief economy (E,D) has at 
least one equilibrium. 
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We will prove the following extension of their result in the next section. 

Lemma 4.1 (Approximation of Games and Markets). For any bimatrix game (A, B), 
let (E,D) = 7r(A,B). Let (u, w) be an e-approximate equilibrium of (E, D) and assume 
u = (x T ,y T ) T and w = (p T ,q T ) T - Then, (x, y) is an O (n^e) -relatively- approximate 
Nash equilibrium for (A,B). 

Lemma 14. 1 1 enables us to prove one of the main results of this paper. 

Theorem 4.2 (Market Approximation is Likely Hard). The problem of finding 
a 1 / 'n ^ 1 ) -approximate equilibrium of a Leontief economy with n goods and n traders is 
PPAD -hard. 

Therefore, if PPAD is not in P, then there is no algorithm for finding an e-equilibrium 
of Leontief economies in time polynomial in n and 1/e. 

Proof. Apply Lemma 14.11 with e = n~ h for a sufficiently large constant h and Corollary 
EH □ 



4.3 The Smoothed Complexity of Market Equilibria 

In the smoothed analysis of the Leontief market exchange problem, we assume that entries 
of the endowment and utility matrices is subject to slight random perturbations. 

Consider an economy with ^E G Kj^pD G R^^Y For a a > 0, a perturbed economy 

is defined by a pair of random matrices (A^A- ) where 6fj and SPj are independent 
random variables of magnitude a. The common two perturbation models are the uniform 
perturbation and Gaussian perturbation. In the uniform perturbation with magnitude a, a 
random variable is chosen uniformly from the interval [— <r, a]. In the Gaussian perturbation 
with variance a 2 , a random variable 5 is chosen with density 




Let D = max (D + A-°, 0) and let E = max (E + A E , 0) . Although we can re- 
normalize E so that the sum of each row is equal to 1, we choose not to do so in favor 
of a simpler presentation. The perturbed game is then given by (E, D). 

Following Spielman and Teng |24| , the smoothed complexity of an algorithm W for the 
Leontief economy is defined as following: Let Tw(E,D) be the complexity of algorithm 
W for solving a market economy defined by (E,D). Then, the smoothed complexity of 
algorithm W under perturbations N a () of magnitude a is 

Smoothediy [n, a] = max ^e^-N^(E) D*-AUD) PV( E > D )] > 

where we use E <— N a (E) to denote that E is a perturbation of E according to N a (E). 

An algorithm W for computing Walrasian equilibria has polynomial smoothed time com- 
plexity if for all < a < 1 and for all positive integer n, there exist positive constants c, k\ 
and &2 such that 

Smoothedvy [n, a] < c • n kl a~ k2 . 



8 



The Leontief exchange economy is in smoothed polynomial time if there exists an algorithm 
W with polynomial smoothed time-complexity for computing a Walrasian equilibrium. 

To relate the complexity of finding an approximate Nash equilibrium of two-person 
games with the smoothed complexity of Leontief economies, we examine the equilibria of 
perturbations of the reduction presented in the last subsection. In the remainder of this 
subsection, we will focus on the smoothed complexity under uniform perturbations with 
magnitude a. One can similarly extend the results to Gaussian perturbation with standard 
deviation a. 

Let (A, B) be a two-person game in which each player has n strategies. Let (E, D) = 
7r(A,B). Let ( A E , A-°) be a pair of perturbation matrices with entries drawn uniformly 
at random from [— a, a]. The perturbed game is then given by E = max (E + A E ,0) and 
D = max (D + A D ,0). 

Let il cr (A, B) be the set of all (E, D) that can be obtained by perturbing 7r(A, B) with 
magnitude a. Note that the off-diagonal entries of E are between and a, while the diagonal 
entries are between 1 — a and 1 + a. 

In the next section, we will prove the following lemma. 

Lemma 4.3 (Approximation of Games and Perturbed Markets). Let (A,B) be a 
bimatrix game with A, B G Rm X 2 ?- For any < a < l/(8n), let (E, D) G II CT (A,B). Let 

(u, w) be an e- approximate equilibrium of (E, D) and assume u = (x T ,y T ) T and w = 
(p T ,q T ) T . Then, (x, y) is an 0{ny/e + n 1 ^^fa)-relatively-approximate Nash equilibrium 
for (A,B). 

We now follow the scheme outlined in [25] and used in [3] to use perturbations as a 
probabilistic polynomial reduction from the approximation problem of two-person games to 
market equilibrium problem over perturbed Leontief economies. 

Lemma 4.4 (Smoothed Leontief and Approximate Nash). Lf the problem of com- 
puting an equilibrium of a Leontief economy is in smoothed polynomial time under uniform 
perturbations, then for any < e' < 1, there exists a randomized algorithm for computing 
an e' -approximate Nash equilibrium in expected time polynomial in n and 1/e'. 

Proof. Suppose W is an algorithm with polynomial smoothed complexity for computing a 
equilibrium of a Leontief economy. Let Tyi/(E,D) be the complexity of algorithm W for 
solving the market problem defined by (E,D). Let N a Q denotes the uniform perturbation 
with magnitude a. Then there exists constants c, k\ and k<i such that for all < a < 1, 



max Ee^jv (S) d^jv (B) [T W (B, D)] < c • n V 



Consider a bimatrix game (A, B) with A,B G ^ x 2 ?. For each (E,D) G iI CT (A,B), by 

Lemma |4.31 by setting e = and a = 0(e'/n 3 ), we can obtain an e'-approximate Nash 
equilibrium of (A,B) in polynomial time from an equilibrium of (E,D). 

Now given the algorithm W with polynomial smoothed time-complexity, we can apply 
the following randomized algorithm with the help of uniform perturbations to find an e- 
approximate Nash equilibrium of game (A, B): 
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Algorithm ApproximateNashFromSmoothedLeontief (A,B) 

1. Let (E,D) = vr(A,B). 

2. Randomly choose a pair of perturbation matrices (A B , A B ) of magnitude 
a. 

3. Let D = max (D + A D , 0) and let E = max (E + A E , 0) . 

4. Apply algorithm W to find an equilibrium (u, w) of (E,D). 

5. Apply Lemma 14.31 to compute an approximate Nash equilibrium (x, y) of 
(A,B). 

The expected time complexity of ApproximateNashFromSmoothedLeontief is bounded 
from above by the smoothed complexity of W when the magnitude perturbations is e'/ra 3 
and hence is at most c • n kl+3k2 (e')~ k2 . □ 

We can use this randomized reduction to prove the second main result of this paper. 

Theorem 4.5 (Hardness of Smoothed Leontief Economies). Unless PPAD C RP, 
the problem of computing an equilibrium of a Leontief economy is not in smoothed polynomial 
time, under uniform or Gaussian perturbations. 

Proof. Setting e' = n~ h for a sufficiently large constant and apply Lemma [4.41 and Corollary 

m □ 



5 The Approximation Analysis 

In this section, we prove Lemma f4.3l Let us first recall all the matrices that will be involved: 
We start with two matrices (A, B) of the bimatrix game. We then obtain the two matrices 
(E, D) = vr(A, B) of the associated Leontief economy, where E = l2 n and 

We then perturb (E,D) to obtain (E,D). We can write D as: 

where for all Vi,j, Zij,mj S [0, a] and aij — <lij,bij — bij £ [—a, a], Note also because 
A, B G ^ji^™ an d < a < 1, A and B are uniform perturbations with magnitude a of A 
and B, respectively. Moreover, Zij and with probability 1/2 and otherwise, they 

are uniformly chosen from [0, o~\. 

Now, let (u, w) be an e-approximate equilibrium of (E,D) and assume u = (x T ,y T ) T 
and w = (p T ,q T ) T , where all vectors are column vectors. 

By the definition of e-approximate market equilibrium, we have: 

Zx + Ay < (1 + e) • 1 

Bx + Ny < (l + e)-l (5) 
(1 - e)E T w < diag(u)D T w < (1 + e)E T w, 
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where diag(u) is the diagonal matrix whose diagonal is u. Since the demand functions are 
homogeneous with respect to the price vector w, we assume without loss of generality that 

ll w lli = IIpIIi + Nli = !• 

We will prove (x, y) is an O (ny/e + n 1 ' 5 y / o : )-relatively-approximate Nash equilibrium 
of the two-person game (A,B). To this end, we first prove the following three properties 
of the approximate equilibrium (u, w). To simplify the presentation of our proofs, we will 
not aim at the best possible constants. Instead, we will make some crude approximations 
in our bounds to ensure the resulting formula are simple enough for readers. 

Property 1 (Approximate Price Symmetry). If ||w|d = 1, < e < 1/2, and < a < 
l/(2n), then 

1 — e — Ana 1 + e 

2 -Ana ~ MW ' Nl1 " 2^W 

Proof. Recall u = (x T ,y T ) T and w = (p T ,q T ) T . By (JHJ and the fact that the diagonal 
entries of E are at least 1 — er, we have 

(l-e)(l-<7) lipid < l T (diag(x)Z T p + diag(x)B T q) = (Zx) T p+ (Bx) T q 
< 3n<7 Hp^ + (1 + e) Nld , 

where last inequality follows from (Bx) < (1 +e)l (obtained from Equation ©), its simple 
consequence X{ < (1 + e)/(l — a) < 3, \/i (because entries of B are between 1 — a and 2 + a), 
and the fact that entries of Z are between and a. 

Applying, ||q|d = ||w|d — ||p|d = 1 — ||p|d to the inequality, we have 

„ „ x + e !+ e 

P i < 7 77 s 7 s < 



Thus, 



(1 - e)(l - a) + (1 + e) - 3na ~ 2 - Ana 
in -, n n 1 — e — Ana 

|q|d = 1 - lipid > 



2 - Ana 

We can similarly prove the other direction. □ 

Property 2 (Approximate Utility Symmetry). If ||w|d = 1 3 < e < 1/2, and 
< a < l/(8n), t/ien 

(l-e)(l-<r)(l-e-4n<r) (1 + e) 2 + na(l + e)(2 - 4na) 

(1 + e )( 2 + 2cT ) - H x lli ' llyld - (± - o)(l - e - Ana) ' 

Proof. By our assumption on the payoff matrices of the two-person games, 1 < aij,bij < 2, 
for all 1 < i,j, < n. Thus, 1 — a < Oy, b{j < 2 + a. By © and the fact the diagonal entries 
of E is at least 1 — a, we have 

(l- e )(l-<T) Pi (l-e)(l-a) Pi (1 - e)(l - a) (2 - 4na) Pi 

*C j , _ i . , i . ^ / _ ... I, I, I, 



<bi|q> + <Zi|p) " (2 + <t) ||q|d + a ||p|d (2 + 2<r)(l + e) 

where the last inequality follows from Property ^ Summing it up, we obtain, 

(1 - e)(l - a){2 - Ana) (1 - e)(l - er)(l - e - Ana) 

I|X|11 - (2 + 2a){l + e) " Pl11 ~ (l + e)(2 + 2a) : 
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where again, we use Property ^ in. the last inequality. 
On the other hand, from (J5J) we have 

(1 + e) ( ei |w) (l + e)(pi + <7) /.,.,. , A +v , + f ^ 

Xj < — — : — ; ; — ; — r < — ; — ; (using w L = 1 and the property of ej 

(b;|q) + (Zi|p) (b;|q) 

(1 + e)( Pi + a) (1 + e)(pj + a) (2 - 4na) 

" (l-aJUqllj " (l-<x)(l-e-4na) ' 1 ' 

Summing it up, we obtain, 

(1 + e)(2 - Ana) na(l + e)(2 - Ana) 

ll x lli — t\ w^j r IIpIIi + 



< 



(1 - <t)(1 - e - Ana) 11111 (1 - a)(l - e - Ana) 
{I + e) 2 + na{l + e){2 - Ana) 



(1 - a)(l - e - Ana) 

We can similarly prove the bound for ||y || 1 . □ 

Property 3 (Utility Upper Bound). Let s = Zx + Ay and t = Bp + Ny. Let X = 

max{e,ncr}. Under the same assumption as Property^ if Sj < (1 + e)(l — a) — \/A, 

(l + e)(2-Ana)(5V\ + a) 
Xi ~ (l-a)(l-e-Ana) ' 

Similarly, if t, L < (1 + e)(l — a) — vA, i/ten 

(l + e)(2-4no-)(5v / A + cj) 



(1 - a)(l - e - Ana) 
Proof. By ©, we have 

(l-e)(l-<r)||w|| 1 < (1 - e)l T E T w < u T D T w = (s|p) + (t|q) 
= ^2 s jPj + (My) + s iPi 

< (l + e)l T E T w- (i + e )(e i |w) + s iPi 

< (1 + e)(l + na) \\w\\ t - (1 + e)(l - a) Pi + s iPi , 

where the inequality immediately after the second equation follows from the last inequality 
of ©. 

Thus, [(l + e)(l-a)-Si]pi < (l + e)(l + ner) - (l + e)(l-ex) < 2e + 3na. Consequently, 
if (1 + e)(l - a) - Si > V\, then p { < (2e + 3na)/V\ < 5\/X By Equation ® 

(l + e)(2-Ana)(pi + a) (1 + e)(2 - Ana)(5V\ + a) 



(1 - a)(l - e - Ana) ~ (1 - a)(l - e - Ana) 
We can similarly establish the bound for y^. □ 
We now use these three properties to prove Lemma [Ql 
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Proof, [of Lemma 14.3] In order to prove that (x, y) is a ^-relatively approximate Nash 
equilibrium for (A,B), it is sufficient to establish: 



x T Ay > (1 — 6) max x T Ay 

y T B T x > (1 - <$) max y T B T x. 

Ily|li=lly|li 

Let s = Zx + Ay. We observe, 

x T Ay = x T (Ay + Zx- Zx+ (A - A)y) = x T s-x T (Zx+ (A - A) y) 

> x T s - llxH, (llZxl^ + ||(A - A)y||J > x T s - a ||x||, (||x||, + ||y||,) 

T 2cr(l + e) 2 + 2no- 2 (l + e)(2 - Ana) 
~ X S (1 -<t)(1 -e-Ana) '^'^ 

= x T s- 0(<7)||x||,, 

where the last inequality follows from Property [21 

Let A = max(e, na). By Property 13 we can estimate the lower bound of (x|s): 



(7) 



(x|s) = ^XjSj> ^ XiSi 

i=l i:Si>(l+e)(l~a)-V\ 



> 



(l + e)(l-<j)- 



> H^i -o[Vx)) [l 

1-0 (nVX + na 



(l +e)(2-4no-)(5y / A + (r) 
(l-<r)(l-e-4ra<r) 

(l + e)(2-4na)(5V\ + a) \ 
(l-<r)(l-e-W)[|x||, I 



x 



On the other hand, by ©, we have Ay < (1 + e)l and hence 



max x 



H x lli=ll x lli 
Therefore, 



Ay = llxH, || Ay|| < (1 + e + a ||yy < (1 + e + 0{a) 



x'Ay > x T s - O(o-) ||x||, 



> Nil 
> 



(1-0 (nVX + na) - 0(a) 



1 



l + e + 0(a) 



1 - 0(nV\ + 



na 



max x Ay 



X , = X 



1 - o(n^/\ + na max x Ay. 

V ll*[|i=l|x[|i 



We can similarly prove 



y T B T x = ( 1 - O (n^/X + na) ) max y T B T x. 

v v )> HyllHIylli 

We then use the inequalities \/A = Y / max(e, ncr) < v / e~+no r < -y/e + y/na and cr < y^cr to 
complete the proof. 

□ 
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6 Remarks and Open Questions 



In this section, we briefly summarize some remarkable algorithmic accomplishments in the 
computation of (approximate) market equilibria obtained prior to this work. We then 
present some open questions motivated by these and our new results. 

6.1 General Market Exchange Problems and Algorithmic Results 

In our paper, we have focused on Leontief market exchange economies. Various other market 
economies have been considered in the literature |231 |§| . 

An instance of a general market exchange economy with n traders and m goods is 
given by the endowment matrix E together with n utilities functions v%, . . . , v n . Then an 
equilibrium price vector is a vector p satisfying 

3 X G M mxn ,Xl < El, X T p < E T p, and Uj(xj) = max{u,(x^) : (x£|p) < (e,-|p)} . 

The pairs (X, p) and ([ui(xi), ...,u n (x n )] ,p) is as also referred to as an equilibrium of the 
exchange market (E, (ui, . . . ,u n )). 

The celebrated theorem of Arrow and Debreu states that if all utility functions are 
concave, then the exchange economy has an equilibrium. Moreover, if these functions are 
strictly concave, then for each equilibrium price vector p, its associated exchange X or 
utilities is unique. 

A popular family of utility functions is the CES (standing for Constant Elasticity of 
Substitution) utility functions. It is specified by an m x n demand matrix D. The utility 
functions are then defined with the help of an additional parameter p £ (— oo, 1] \ {0}: 



As p —* — oo, CES utilities become the Leontief utilities. When p = 1, the utility 
functions are linear functions. 

Remarkablely, an (approximate) equilibrium of an exchange economy with linear utilities 
functions can be found in polynomial time |20| 1121 127| I13j . In fact, Ye shows that an e- 
equilibrium of such an market with n traders and n goods can be found in 0(n 4 log(l/e)) 
time '27 ]■ If data is given as rational numbers of L-bits, then an exact equilibrium can be 
found in 0(n 4 L) time. 

A closely related market exchange model is the Fisher's model |23| . In this model, 
there are two types of traders: producers and consumers. Each consumer comes to the 
market with a budget and a utility function. Each producer comes to the market with an 
endowment of goods and will sell them to the consumers for money. An equilibrium is a 
price vector p for goods so that if each consumer spends all her budget to maximize her 
utilities, then the market clears, i.e., at the end of the exchange, all producers sold out. 



Even more remarkablely, an (approximate) equilibrium in a Fisher's economy with any 
CES utilities can be found in polynomial time [TTI 1271 l28l fTTH IT4] . 
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6.2 Open Questions 

Our results as well as the combination of Codenotti, Saberi, Varadarajan, and Ye [H] and 
Chen and Deng [3] demonstrate that exchange economies with Leontief utility functions are 
fundamentally different from economies with linear utility functions. In Leontief economies, 
not only finding an exact equilibrium is likely hard, but finding an approximate equilibrium 
is just as hard. 

Although, we prove that the computation of an O(l/n 0< - 1 ^)-approximate equilibrium of 
Leontief economies is PPAD-hard. our hardness result does not cover the case when e is a 
constant between and 1. The following are two optimistic conjectures. 

Conjecture 4 (PTAS Approximate LEONTIEF). There is an algorithm to find an 
e-approximate equilibrium of a Leontief economy in time 0(n k+e c ) for some positive con- 
stants c and k. 

Conjecture 5 (Smoothed LEONTIEF: Constant Perturbations). There is an algo- 
rithm to find an equilibrium of a Leontief economy with smoothed time complexity 0(n k+cr c ) 
under perturbations with magnitude a, for some positive constants c and k. 
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